Phase dependence of localization in the driven two-level model 
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A two-level system subjected to a high-frequency driving field can exhibit an effect termed "coher- 
ent destruction of tunneling" , in which the tunneling of the system is suppressed at certain values 
of the frequency and strength of the field. This suppression becomes less effective as the frequency 
of the driving field is reduced, and we show here how the detailed form of its fall-off depends on the 
phase of the driving, which for certain values can produce small local maxima (or revivals) in the 
overall decay. By considering a squarewave driving field, which has the advantage of being analyt- 
ically tractable, we show how this surprising behavior can be interpreted geometrically in terms of 
orbits on the Bloch sphere. These results are of general applicability to more commonly used fields, 
such as sinusoidal driving, which display a similar phenomenology. 
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Since the pioneering work of Hanggi and co-workers* 
it has been known that a quantum system subjected to a 
periodic driving field can exhibit an unusual effect termed 
"coherent destruction of tunneling" (CDT), in which the 
tunneling dynamics of the system are quenched at vari- 
ous specific values of the strength and frequency of the 
applied field. One of the simplest systems in which this 
effect appears is the driven two-level system*, which has 
been applied to describe many physical situations, such 
as an electron moving in a double- well potential*, molec- 
ular switches*, the single-Cooper-pair boji* and super- 
conducting tunnel junctions 6 . Controlling the dynam- 
ics of such systems has become of increasing importance 
due to possible technological applications in quantum 
computing, and the CDT effect provides one method of 
achieving such control which preserves the quantum co- 
herence of the system. 

The periodicity of the driving field allows us to invoke 
the Floquet theorem, and expand solutions of the time- 
dependent Schrodinger equation* as 
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where tj is termed the quasienergy, and 4>j (t) is a function 
with the same periodicity as the driving field, called a 
Floquet state. A necessary condition for CDT to occur is 
that the quasienergies must be either degenerate or close 
to degeneracy*. If we neglect the time-dependence of the 
Floquet states it is clear that the dynamics of the sys- 
tem will indeed be frozen when this condition is satisfied. 
Neglecting the intrinsic time-dependence of the Floquet 
states is a reasonable approximation when the frequency 
of the driving field is high (the "high-frequency limit"—), 
as the period of the field is then much shorter than the 
typical timescale for tunneling processes, given by the 
Rabi period. At lower frequencies, however, the Floquet 
states can have a non-trivial time-dependence, and the 
degree of CDT can be markedly reduced if the Floquet 
states themselves have a large amplitude of oscillation 8 . 
In this work we will show that fall-off in CDT produced 
by this effect has a surprising dependence on the phase 
of the driving field. Although phase-effects have been 



briefly noted in previous work*, we provide here the first 
systematic study and explanation of this effect. We shall 
show how the phase can be chosen to maximize the degree 
of CDT, and also how certain choices of phase can pro- 
duce a non-monotonic decay of CDT, containing small 
local maxima or "revivals" , as was seen previously in 
Ref.ia. 

We consider a two-level system described by the Hamil- 
tonian 
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where <7, are the standard Pauli matrices. We parame- 
terize the driving field as Ef(t), where E is the strength 
of the driving, and f(t) is a T-periodic function of unit 
amplitude and frequency u>, which describes the wave- 
form of the driving field. The eigenstates of the undriven 
system consist of two extended states: namely a symmet- 
ric ground state separated by the level splitting, A, from 
the anti-symmetric excited state. To investigate the tun- 
neling dynamics of the system it will prove useful to also 
define localized states formed by the sum and difference 
of these eigenstates, which we label as \L) and \R). The 
degree of tunneling can be conveniently assessed by ini- 
tializing the system in state \L) and then evolving it in 
time under the influence of the Hamiltonian (J2J) . By mea- 
suring the probability that the system remains in state 
\L), Ph(t) — \(4>(t)\L)\ 2 , we can quantify the degree to 
which tunneling is suppressed by finding the minimum 
value of Pl it) attained, which we term the "localization" . 
When the tunneling is completely destroyed Pt(t) does 
not change with time, and thus the localization is equal 
to one. Conversely, if the tunneling is not completely 
suppressed, the localization will take a lower value, and 
will be zero if the particle is able to completely tunnel 
from \L)to\R). 

We show in FigQ] the localization produced as the 
frequency of the driving field is reduced, obtained by 
the numerical evolution of the system using a Runge- 
Kutta techniquei*, for two types of driving field: a sine 
wave and a squarewave. For both fields the driving 
strength was chosen so that the quasienergies of the 
system were degenerate and so CDT occurs. As was 
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FIG. 1: Localization produced by (a) a sinusoidal driving 
field and (b) a squarewave driving field. Solid line indicates 
the result for a driving with zero phase-shift, the dotted line 
for a phase-shift of 7r/4, and the dashed line for a phase-shift 
of 7r/2. For high-values of ui the localization is excellent in all 
cases, but drops rapidly for A/lo > 0.5. For the phase-shift 
of 7r/2 the localization goes through a small local minimum 
before reaching zero. Insets: A single period of the zero phase- 
shifted driving fields f(t). 



shown in Ref. n », this degeneracy condition defines a set 
of smooth one-dimensional curves in (E, u>) space called 
"crossing-manifolds" 3 . When the quasienergies are de- 
generate it is clear from Eq. (JIJ that the time-dependence 
of the system is given solely by the Floquet states. Since 
these states are periodic, it is sufficient to evolve the sys- 
tem over just one period of the driving to study its proper- 
ties, giving a considerable saving of computational effort. 

Although the two driving fields are very different, they 
produce strikingly similar forms of localization. In both 
cases in the high-frequency regime (A/cj <C 1) the local- 
ization has a value of one, irrespective of the phase of the 
driving field. This indicates that the tunneling is com- 
pletely destroyed and thus the system remains frozen in 
its initial state. As the value of u> is reduced, the par- 
ticle is able to tunnel more freely, and the degree of lo- 
calization falls. It can also be clearly seen that at low 
frequencies the phase of the driving field has a signifi- 
cant effect on the form of this fall-off. When the field 
has no phase-shift (shown in the insets to FigQJ the lo- 
calization decays monotonically to zero as the frequency 
is reduced. The initial effect of increasing the phase of 
the driving field from zero is to cause the localization to 
drop more rapidly. For example, at A/lo = 1 the lo- 
calization produced by a ir/4 phase-shifted field is half 
that of the zero phase-shift case, while for a phase shift 
of 7r/2 the localization is about four times smaller. A 
further effect occurs for values of the phase shift close to 
its maximum value of ir/2, in which after rapidly decay- 
ing to zero, the localization subsequently passes through a 
small local maximum, before again vanishing at the same 



frequency as the zero phase-shift case. 

To understand how this behavior comes about, it is 
useful to visualize the dynamics of the two-level system 
geometrically, by making use of the Bloch sphere rep- 
resentation well-known from quantum optics and NMR 
studies, and now commonly used in quantum informa- 
tion theorjii^. Any pure state of the system can be rep- 
resented by a point on the surface of the Bloch sphere, 
and we can note in particular that the localized states 
\L) and \R) correspond to the points [0,0, ±1]. We 
shall consider just the case of squarewave driving, as 
it displays the same phase-dependent behavior as sinu- 
soidal driving, but also has the useful property of al- 
lowing exact solutions for the time-development of the 
system to be obtained. For this form of driving, the 
time-dependent Hamiltonian J2J decomposes into piece- 
wise constant parts, H± — (A/2)a x ± (E/2)a z . These 
operators can be interpreted conveniently as an interac- 
tion between the Bloch vector s and a fictitious magnetic 
field directed along the axes r± = [A,0,±E], and thus 
the time-evolution of the system can be viewed as suc- 
cessive Larmor rotations about these axes. This insight 
is the key to a simple method for interpreting and under- 
standing the dynamics of the system. 

The operators H± can be exponentiated straightfor- 
wardly, yielding the following result for the unitary op- 
erator describing the evolution of the system under the 
influence of H±: 



u±(t) = /cos nt- 



sin fit . . „ 

[A cr x ± E * z ] . (3) 



In this expression fl = \JE 2 + A 2 /2 is the Rabi fre- 
quency. The general time-evolution operator for a 
squarewave driving field can now be expressed simply 
as the product of successive factors of U±. Examination 
of Eq[2| reveals that the propagator for one period of the 
driving, U(T) = U+(T/2)U-(T/2), is equal to the iden- 
tity (and thus the quasienergies are both zero and CDT 
occurs) if sin SIT/ 2 = 0. This indicates that CDT only 
arises if the driving frequency is in resonance with the 
Rabi frequency, f2 = nto where n is an integer. This 
resonance condition may be expressed alternatively as 



(2nY 



(4) 



It is interesting to note that this equation for the crossing- 
manifolds for a squarewave driving field, which describes 
concentric circles in (E/uj,A/uj) space, is identical to 
that hypothesized from purely numerical evidence in 
Ref.il. 

We now consider the explicit time evolution of the sys- 
tem, and begin with the simplest case of zero phase- 
shift. The state \L) in which the system is initial- 
ized corresponds to the north-pole of the Bloch sphere 
(s = [0,0,1]). During the first half-period of the driv- 
ing, this initial state evolves under the influence of H + , 
making a Larmor rotation about the axis r + . The res- 
onance condition between the driving frequency to and 



A/0) 



FIG. 2: Orbits on the Bloch sphere for driving parameters 
A/lo = 1/2, E/uj = a/T5/2, i.e. on the first crossing-manifold 
(n = 1). The left plot shows the time-evolution for zero phase- 
shift, and on the right for a phase-shift of 7r/2. 



the Larmor frequency means that during this interval 
the Bloch vector makes n complete revolutions, and thus 
at t = T/2 the system returns to its initial state. Dur- 
ing the second half-period a similar process occurs, with 
the Bloch vector making n revolutions about the axis r_ . 
The system thus traces out a symmetric "figure-of-eight" 
on the surface of the Bloch sphere, as shown in Fig|21 In 
the high-frequency limit, when w» A and thus E ~ lu 
from Eq0| the rotation axes are almost parallel to the 
z-axis, and so the loops traced by the Bloch vector have 
a very small radius. As a result, the orbit remains close 
to \L) at all times, and the system exhibits a high de- 
gree of localization. As E is reduced, however, the angle 
between the rotation axes increases and the loops traced 
by the Bloch vector become larger. The system can thus 
explore further around the Bloch sphere and approach 
closer to \R), thereby reducing the localization. It is a 
simple matter to explicitly derive an expression for Pl (t), 
from which the following result for the localization can 
be obtained 



Pj am = PL(T/4) = l- 



4n 2 



(5) 



In FigEJit can be seen that this result indeed correctly de- 
scribes the fall-off in localization for a zero phase-shifted 
field. 

We now consider the effect of the phase-shift. For a 
phase-shift of tt/2 the system first evolves under H + for 
a quarter-period, then under if_ for an interval of T/2, 
and completes the final quarter-period under H+. Ac- 
cordingly, during the first quarter-period of the driving 
the Bloch vector makes n/2 revolutions about r + , and it 
is thus now important to distinguish whether n is even 
or odd. In the former case, the Bloch vector will make 
an integer number of revolutions about r + in the first 
quarter-period, and so it will end this interval at \L). It 
will then make n complete revolutions about r_ during 
the next half-period, and will finally complete its motion 
by making n/2 revolutions about r+. For even values of n 
the Bloch vector thus traces out the same figure-of-eight 
orbit as for the zero phase-shift case, and so although the 
order in which it traverses the loops is different, identical 



FIG. 3: Comparison of the localization on the first crossing- 
manifold (n = 1) with the analytical expressions derived in 
the text. EqQ3is denoted by the solid line, Eq^Jby the dashed 
line. It can be seen that for zero-phase shifted driving (hollow 
circles) the localization is described by EqE] For 7r/2 phase- 
shifted driving (solid circles) the localization is given by min 
[Pl(T/4),Pl(T/2)], and the crossover between these quanti- 
ties at A/uj — V3 produces the cusp- like revival feature. 



values of localization will be produced. 

If n is odd, however, a different behavior will occur. 
At the end of the first quarter-period the system will not 
have returned to \L), as the Bloch vector will only have 
made a half-integer number of rotations about r + . The 
Bloch vector has thus been rotated away from the r_ axis, 
and so under the influence of H_ the Larmor rotation 
it makes around this axis will have a different radius, as 
shown in FigEl The second loop clearly approaches much 
closer to \R) than for the case of the figure-of-eight orbit, 
and thus the localization is reduced. 

In Fig^we show this effect in detail for driving param- 
eters on the first-crossing manifold (n = 1), by displaying 
the Bloch sphere orbits projected onto the (s y , s z ) plane. 
For A/u> = 1 the system makes a similar orbit to that 
shown in Fig [3 The minimum value of localization (i.e. 
the closest approach to \R)) occurs at t = T/2, in con- 
trast to the case of zero-phase shift, for which the min- 
ima occur at t = T/4 and 3T/4. Although it is possible 
to calculate explicit expressions for Ph(t), these are in 
general somewhat cumbersome. We instead present for- 
mulae just for times that are multiples of T/4, which, as 
can be seen from Fig^J correspond to the turning points 
in P L (t). 

During the first and last quarter-periods of driving the 
time-evolution of the system is identical to that of the 
case with zero phase-shift, and accordingly Pl(T/4) is 
given by Eq|5] It is also straightforward to show that 



P L (T/2) 



(A/u) 
2n 2 



■ sin 2 nn/2 



(6) 



In Fig |3| it can be seen that this expression indeed cor- 
rectly describes the rapid fall of localization for the tt/2 
phase-shifted field, which reaches zero at A/uj = \[2. At 
this value of frequency, the two rotation axes make an 
angle of 7r/4 to the z-axis, producing the highly symmet- 
ric time-evolution shown in Fig^J}, in which the system 
passes through \R) at t = T/2. As A/w is increased fur- 
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FIG. 4: Time-evolution of the system in the low-frequency 
regime u < A. On the left, the evolution on the Bloch sphere 
projected onto the (s y , s z ) plane. The solid line indicates the 
evolution of the system over the first half-period, the dotted 
line over the second half-period. The upper (lower) crosses 
indicate the states \L) (|-R)). On the right, the time-evolution 
of the probability PL(t) over one period of the driving. The 
solid line indicates zero phase-shifted driving, the dashed line 
has a phase-shift of n/2. 



ther, however, the trajectory of the orbit no longer passes 
exactly through \R) and thus the localization increases, 
producing the "revival" effect. This increase continues 
until A/lj = V3, at which P L (T/4) and P L (T/2) become 
equal. As shown in Fig0Ji, at this value of frequency the 
system does not display any time evolution during the in- 
terval (T/4, 3T/4), since the action of first quarter-period 
of the driving is to rotate the initial state to align with 
r_ . The following half-period of driving under _ff_ , thus 
does not alter the state of the system. 

For still larger values of A/w, the rotation about the r_ 



axis acts to move the system away from \R) (see FigQJl), 
and consequently the minimum value of Pl no longer 
occurs at T/2, but at T/4 instead. In this regime the 
localization produced by the 7r/2 phase-shifted field is 
therefore identical to that produced by the zero phase- 
shifted case. It can be clearly seen in FigOHthat it is this 
crossover in behavior that gives rise to the revival feature 
seen in Fig 2] 

In summary, we have studied how the CDT effect de- 
pends upon the frequency of the driving field. In the 
high-frequency regime, the positions of the quasienergy 
degeneracy can be found analytically 2 ^ and the degree 
of CDT is not affected by the phase of the driving field. 
As the frequency is reduced a more complicated picture 
appears: the quasienergies drift away from their high- 
frequency values along one-dimensional manifolds, and 
the degree of CDT induced depends not only on the val- 
ues of field strength and frequency, but also on its phase. 
We have firstly shown that the form of the crossing- 
manifolds for the squarewave driving field (Eq0J that 
was conjectured previously on numerical evidenced is in 
fact exact. By using the Bloch sphere representation we 
have given a simple method of understanding how the lo- 
calization falls as the frequency of the driving is reduced, 
and have analytically confirmed the empirical observa- 
tion in Ref£ that the localization is maximized for zero 
phase-shifted driving. We have also clarified how the re- 
vival features seen for it/ 2 phase-shifted driving^ arise, 
and have demonstrated that for crossing-manifolds with 
even values of n that this phenomenon does not occur. 
This opens up new prospects for experiment in the low- 
frequency regime, which is both easier to attain and also 
has the advantage that the driving field is less likely to 
drive transitions to higher energy-states (thereby break- 
ing the two- level approximation). 
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